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ABSTRACT 

We study the treatment of the constraints in stochastic quantization method. We 
improve the treatment of the stochastic consistency condition proposed by 
Namiki et al. by suitably taking account of the Ito calculus. Then we obtain an 
improved Langevin equation and the Fokker-Planck equation which naturally leads 
to the correct path integral quantization of the constrained system as the stochas- 
tic equilibrium state. This treatment is applied to 0{N) non-linear a model and 
it is shown that singular terms appearing in the improved Langevin equation can- 
cel out the S^{0) divergences in one loop order. We also ascertain that the above 
Langevin equation, rewritten in terms of independent variablesis, actually equivalent 
to the one in the general-coordinate-transformation-covariant and vielbein-rotation- 
invariant formalism. 
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1. 1. Introduction 



Stochastic quantization method (SQM) was first proposed by Parisi and Wu. They| 
showed that the method could be apphed to gauge theory without the gauge fixing 
procedure. That is, in SQM it is not necessary to introduce the Faddeev-Popov 
ghost fields. Nevertheless the method produces the same contributions as those due 

[3] 

to ghost fields, which was perturbatively confirmed first for Yang-Mills field and 
recently for non-Abelian anti-symmetric tensor field!^' 

How to handle the constrained system in SQM was discussed by Namiki et al. 
in ref.5. They constructed Langevin equation for the system under the holonomic 
constraints by imposing the stochastic consistency condition. In the path integral 
quantization method the constraints introduce a determinant factor into path in- 
tegral measure, which reqinres that in SQM Langevin equation for the constrained 
system is constructed so that the equilibrium Fokker-Planck distribution derived 
from the Langevin equation has the same determinant factor. They showed that 
the equilibrium distribution coincided with the path integral distribution. Neverthe- 
less they had to use the 5-dimensional stochastic path integral representation of the 
transition probability distribution of stochastic process and could not derive Fokker- 
Planck equation directly from Langevin equation, because they did not take account 
of Ito calculus in their treatment of stochastic consistency condition. One of the 
main purposes of this paper is to improve Langevin equation of ref.5 by suitably tak- 
ing accoimt of Ito calculus, and to show that the improved Langevin equation leads 
to the Fokker-Planck equation which directly gives the correct path integral repre- 
sentation as the equilibrium distribution. Our improvement introduces some singular 
terms proportional to S{Q) in the Langevin equation. So, the Langevin equation of 
ref.5 is correct in the dimensional regularization scheme, but not in other regulariza- 
tion schemes and the general path integral distribution with constraints cannot be 
directly obtained as the equilibrium Fokker-Planck distribution. 

On the other hand, if the constraints are solved explicitly and the system is de- 
scribed in terms of the independent variables only, the action of the system proves 
to have generally field-dependent metric. In this case we must apply the general- 
coordinate-transformation (GCT)-covariant and vielbein-rotation (VR)-invariant 
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Langevin equation to the system. It was not clear whether the latter Langevin equa- 
tion is equivalent to the above Langevin equation for the constrained system. This 

point will be clarified in this paper. 

Next we discuss about 0(N) non- linear a model as an example of such con- 
strained system. The model was studied by many authors.'^^'We apply the above 
two methods to the model and clarify the role of singular terms introduced into 
the Langevin equation with the stochastic perturbation theory. It will be shown 
that in both these methods these singular terms are necessary to cancell the 5"^(0) 
divergences appearing in one-loop expansion. 

This paper is organized as follows. In section 2 we study the improved treatment 
of the constraints in SQM and show directly that the equilibrium Fokker-Planck 
distribution coincides with the path integral distribution. In section 3 it is shown 
that the improved Langevin equation for the constrained system is equivalent to the 
GCT-covariant and VR-invariant Langevin equation. In section 4 we apply both 
the improved Langevin equation for the constrained system and the GCT-covariant 
and VR-invariant Langevin equation to 0(N) non- linear a model and examine the 
cancellation of (5"'(0) divergences. In section 5 we give conclusion and summary. 
In Appendix, we ascertain that an assumption, which is introduced in section 3, is 
satisfied at least in 0(N) non-linear a model. 

2. 2. Constrained system in SQM 

In this paper we consider the system with variables = 1, 2, • • • , A^^), regular 

Lagrangian (/^ = 1, 2, ■ ■ ■ , n) and a set of constraints 



In the path integral quantization method the transition amplitude is given by 



Fa{qi) = 0, 



(a=l,2,... 



M;N > M). 



(2.1) 




DqDJ ^JdetDab exp[- / rx{L{q,dq) - JaFa)]. 



(2.2) 



Dab 



dFg dFb 
dqi dqi 



(2.3) 
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Following the method "time by time constraint" proposed in ref.5, we quantize 
the above singular system in SQM. The treatment of the stochastic consistency 
condition is improved by taking account of Ito calculus. It will be shown that the 

improvement is essential for the acquisition of the Fokker-Planck equation which 
directly leads to the path integral representation as the equilibrium state. 

For the above system, Langevin equation is 

dS OF 

dqi{t) = qi{t + dt) - qi{t) = -—dt - ^Kdt + dWi, (2.4) 



S = J d'^xL, (2.5) 
where Xa is Lagrange multiplier and dWi is defined as 

dWi{t) = Wi{t + dt) - Wi{t), (2.6) 

t 

Wi{t) = J dt'riiit'), (2.7) 

{dWi{x, t)dWj{x', t)) = 26ijS{x - x')dt, (2.8) 

(r^,(x, t)r],{x', t')) = 2SijS''ix - x')Sit - t'). (2.9) 

and called Wiener process. From (2.8) we may regard dW as order \/dt. Lagrange 

[5] 

multiplier is determined by the stochastic consistency condition 

Fa{q{t)) = Q, (2.10) 
where dot denotes fictitious-time derivative. Besides, we demand the initial condition 

Fa{q{t = to)) = (2.11) 
in order to have F = at any t. In SQM, for fictitious-time derivative of any function 
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f{Q) 



[6] [11] 



pm) + 1^, (2.12) 

[12] 

unless we use Stratonovich calculus. In ref.5 the above fact is not taken into account. 
(Their Langevin equation is correct only in the dimensional regularization scheme. 

In other regularization schemes, the Langevin equation is incorrect and the Fokker- 
Planck distribution does not coincide with the path integral distribution.) If \dt in 
eq.(2.4) does not contain terms of order \fdt^ dq = dW to order ^/dt. From (2.4) we 
get 

= t-<%+^" ''''+''^'>+2(i-2'')s^^''»*''^^^^+o«^)')- (2-13) 

rgi 

to order dt. Here we use the generalized Ito formula and the product of Wiener 
process dW and any function f{q) is defined as 



{f{q) dW}{t) ^ {bf{q{t + dt)) + (1 - b)f{q{t))}dW{t), 



Q<b<l, (2.14) 

where b = ^ corresponds to Stratonovich calculus and & = to Ito calculus. Requir- 
ing dFa = 0, we get 

= D-J{^i-s^dt + m) + 2(1- ^)g;^dWm,}, (2.15) 



Here Xadt contains D~^^^dWi of order which is inconsistent with the above 
assumption that Xadt does not contain terms of order ^/di. So we assume alterna- 
tively that the Xadt term contains terms of order Vdi like eq.(2.15). Then, from 



-5 



(2.4) 

dF 

dqi = dWi - —^Xadt 
to order \/di. From (2.17) the modified expression of dF is 



Kij = Sij - Rij, Rij = 'o^^ab'o^^ (2-18) 

to order dt. Here Rij, Kij are projection operators, vertical each other. The final 
term in RHS of (2.13') does not exist in ref.5. From the consistency condition (2.10) 
the correct expression of Xdt becomes 

SQ 1 B^F 

^""'^ ^ ^-^'"^ + "^^'^ + 2 " ^^^d^j^^'^'^^'^^^''^^'^- ^^-^^'^ 

The above expression is surely correct, because, with the help of (2.15'), dq has the 
same terms as (2.17) to order ^/di. From (2.8) and (2.15'), Langevin equation (2.4) 
becomes 

If we multiply (2.19) by we obtain the expression (2.13') or (2.10). Therefore, 
the consistency condition (2.10) is embedded in eq.(2.19) and eq.(2.19) means N- 

71^ independent difi:erential equations. In ref.5 Langevin equation did not have the 
singular term in (2.19). The same equilibrium Fokker-Planck distribution as (2.2) 
could not be derived directly from the Langevin equation of ref.5, while, due to 
the singular term, we can derive the correct Fokker-Planck equation directly from 
Langevin equation (2.19) as shown below. 
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In order to construct the Fokker-Planck equation, we introduce the expectation 
value of fictitious-time derivative of arbitrary function g{q) 

{g{qm= J Dq g{q)P{q,t), (2.20) 

where P{q, t) is the transition probabihty distribution. Using integral by parts and 
generalized Ito formula, we obtain the Fokker-Planck equation 

due to the singular term in (2.19). Eq.(2.21) cannot be derived from the Langevin 
equation of ref.5. The probability distribution can include any function f{Fa) be- 
cause ^ is vertical to the projection operater Kij. In the equilibrium limit t — > oo, 
the probability distibution must satisfy P = 0. In the limit the equation has a 
solution 

P{q) = J DvaVdeWabexp{-S- J dl'x VaFa), (2.22) 

where we chose J Dvexp{—VaFa) as f{Fa) in accordance with the initial condition 
(2.11) in Minkowski space. The above equilibrium probability distribution coincides 
with eq.(2.2). However, it is strange that the Langevin equation has divergent term 
of 5"(0). In section 4 we examine 0(N) non- linear a model as an example of the 
system under constraint (2.1). In the same example we also show perturbatively 
that the term proportional to 6'^{0) in eq.(2.19) is needed. 

3. 3. Equivalence to GCT-covariant 
and vr-invariant langevin equation 

In this section we shall show that the improved Langevin equation (2.19) is 
actually equivalent to the Langevin equation in the GCT-covariant and VR-invariant 
formalism where the constraint (2.1) is explicitly solved and the equation is expressed 
in terms of independent variables only. 
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In general the system may have field-dependent metric (or kernel) Gab{q) after 
the constraint is solved. According to ref.7, a system with field-dependent metric is 
described by the GCT-covariant and VR-invariant Langevin equation 



= X^dt + E^dW"^, (3.1) 



where 



= A,5,m,n = (1,2,- ■■,iV-M; x). (3.2) 

G^^ is the inverse of metric Gab, the summation with respect to B includes space- 
time integration and dW^ is Wiener process defined in section 2. GCT-covariance 
and VR-invariance mean that Langevin equation is transformed covariantly under 
general coordinate transformation q ^ q' — f{q) and is invariant under vielbein 
rotation E"^ E\k^^. dq^ and E^dW^ are not transformed covariantly and 
two extra terms in (3.2) are required to be GCT-covariant and VR-invariant. 

In order to decompose variables into constraint variables and independent ones, 
we introduce a new set of variables {Q^}{ijl = 1, 2, • • • , A^).'^'(5^'s are expressed in 
terms of qis {i = 1, 2, • • • , N) as 

5Q^ = e^,5g,, or ^ = e% (3.3) 

where e^- is vielbein field defined as follows. First e"^ and ea,i (a = N—M+1, ■ ■ ■ ,N) 
are defined as 

^~ ^q^' "^"'^ " dqi ' ^^^^ 

i.e. = Fa- Then, e"^i {A= 1,2, ■ ■ ■ , N — M) is chosen so as to satisfy 



e^ea,i = 0, (3.5) 
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and its inverse e^^i is defined as 



eA,r - (^~')Aijei = eiei (3.6) 

where we assume that g^^ is non-singular. From the above definition it turns out 
that e'^^ and e^^^ satisfy the following relations 

ele''^ = e^ea,i = eA,ie''i = e^.iea.i = 0, (3.8) 

Kije'^j = KijCaj = 0, Rije"] = RijCAj = 0. (3.10) 

dei(e^.) 7^ 0. (3.11) 

From (3.11) the manifold spanned by g^'s is identical with the one by Q^'s. With 
the help of the same discussion as made about eq.(2.13'), dQ^^ is written as follows 



dQ^' = Q^(i + dt) - Q^(i) = -^dqi + (1 - 2b)—^K,jdtS''{Q), 

= e^'idqi + (1 - 2b)-^e^ieA,jdtS''{0). 

Then, from (2.13') constraint variables Q^'s satisfy Langevin equation 



(3.12) 



+ {l-2b)^eieA,jdtS-{0), ^^"^^^ 

= 0, 

to order dt. From (3.13) and the initial condition (2.11), constraint variables (5"'s 
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are zero for all t. As for the independent variables Q^^s we get 

+ (l-26)-^e%,,di^"(0), (3.14) 
oqj 



The above Langevin equation is not invariant under vielbein rotation e"^ — ^ e^-Pv' ^{Q) 
because Wiener process dW^ is defined in a manifold spanned by original variables 
gi's and we must not consider the rotation in the manifold. If we perform field- 
dependent rotation in the manifold, -dW^ is not Wiener process, i.e.(A*^dM^-^) ^ 0. 
In order to reduce (3.14) to the form of (3.1) we decompose the vielbein as follows: 

A = E'\{Q)c\{Q), 7=(l,2,...,Ar-M), (3.15) 



e',e\ = 5'\ eM = i5W'=/^ (3-16) 

and define dW^ by 

dW' ^ e\dW^ - 2h^E^je'r{^)- (3-17) 

Then, we obtain 

{dW^) = 0, {dW\x)dW\y)) = 28^'^dt5'^{x - y). (3.18) 
dW^ is desirable Wiener process and with dW^ Langevin equation is written as 

dQ^ = -/5_dt+_(E46^,)E56^,di5-(0)-26^EV^dic5-(0)+^4rfW^'. 

(3.19) 

From now on we regard E"^j as vielbein corresponding to E^^^ in (3.1). In fact the 
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above Langevin equation is invariant under vielbein rotation 



E^j E^jA-^i. (3.20) 

Langevin equation (3.19) is a little different from GCT-covariant and VR- invariant 
Langevin equation (3.1). Here we assume 

where is covariant derivative in Riemanian manifold spanned by Q^. Eq.(3.21) 
is usually presumed, because g^^ satisfies the metric condition VaS^^ = and the 
latter condition makes eq.(3.21) naturally understandable. We ascertain in Appendix 
that the above assumption is satisfied in 0(N) non-linear a model. With (3.21), 
Langevin equation (3.19) is reduced to 

= -/''^^*+;)^a^(v^/'')^"(0)^*-2&||i^^V'c5"(0)di^ 

(3.22) 

g = det{gAB). (3.23) 

Covariant derivative in (3.21) does not include spin connection because, as mentioned 
above, the rotation must not be considered in the manifold spanned by qis. Eq.(3.22) 
is equivalent to GCT-covariant and VR- invariant Langevin equation (3.1). Thus it 
is ascertained that Langevin equation (2.19) for the constrained system is equiva- 
lent to GCT-covariant and VR- invariant Langevin equation (3.1). That Langevin 
equation also has divergent term including 5"(0) and, in the next section, we apply 
the Langevin equation to 0(N) non- linear a model and see the cancellation of 5"(0) 
divergences. 
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4. 4. 0(N) NON-LINEAR (7 MODEL 

0(N) non- linear u model is defined by action 



S=\j (Pxd^^id^^i, (//=l,2,---,n;i = l,2,---,7V), (4.1) 



and constraint 



F = ^i^i - - = 0, (4.2) 
a 



where a is constant. The constraint is an example of (2.1). Applying Langevin 
equation (2.19) to the model, we obtain 

6i = {5ij - a^i^j){d''^j + r]j) - (1 - 2h)a{N - 1)$,5'^(0). (4.3) 

The above equation includes (2.19)-type constraint (4.2) and means N-1 indepen- 
dent equations. From eq.(4.2) $j has non- zero vacuum expectation value. We shift 
the field 

0^ = ^i- < >, (4.4) 
1 



Vi = 



< $i >, ViVi = 1, (4.5) 



and with the shifted field (j) Langevin equation is written 

4)^ = {Kij - y/a{vi(pj + Vjcpi) - a(f)i(f)j}{d'^(pj + r]j} 



■(1 - 2b)a{N - l){^v^ + 0.)5"(O), (4.6) 



Kij =Sij - ViVj. (4.7) 
Going to momentum space and integrating eq.(4.6) with respect to t, we get 



(f)i{k, t) = j drGijik, t-T)[ Kji r]i{k, r) + Ij{k, r) + Jj{k, r) 



-(1 - 2b)a{N - 1){^ + ct>j{K r)}6-m. (4-8) 



Gij{k, t) = [exp{-K kH}]ij = exp{-kH)Kij + {5ij - Kij), (4.9) 



-j^-y^5''{k - p- q){vi5im + Vmhl)(l>l{p, t){-q^(i>miq, t)+r]miq, t)}, 

(4.10) 

(27r)2^ ^"(^ - P - g - 0^^(P, + Vjir, t)}, 

(4.11) 

After solving the above equation by iteration, we express the result graphically 



Fig. 1 



where we denote rj hy a cross or an encircled cross, G a line and a6{0) a bullet, 
respectively. We calculated the one-loop corrections of the two-point function and 
obtained six 5(0)-divergent diagrams. 



Fig. 2 



Each of the six diagrams contributes respectively 



(2a) 2aNS^{k + k')K,,j^^^^ j dy, 



(26) -4«iV5"(fc + k')e{0)K^jj^^-^-^ j dy, 
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(2c) 2av'S-{k + m^. (fc2 + fc/2)2 / 



(2d) -4a^2^-(^ ^ ;,');^. J dy, 



(2e) 4«^;2c5"(fc + fcO^(0)^»i + fe/2)2 / 

(2/) -2a(l - 2&)(Ar - l)(27r)"5"(0)5"(fc + k')K,j-^j^^^. 

The contributions from all the 5(0)-divergent diagrams cancel out'ii we put ^(0) = 
h. It turns out that 5"(0)-divergent term of eq.(2.19) or eq.(4.3) is necessary to the 
cancellation of the (5(0) divergences in one-loop order. 

Now we apply Langevin equation (3.1) to the model. As discussed in section 3, 
GCT-covariant and VR- invariant Langevin equation is equivalent to eq.(2.19). We 
show that the (5"'(0) divergences really cancel out. 

With substituted by means of constraint (4.2), action (4.1) becomes 

Sn-1 = j rf"xy"^($)<9^$«a^$;3, (4.12) 



g'^^m = + T^^^, a, 7 = 1, 2, • • • , iV - 1. (4.13) 

We can use Langevin equation (3.1) for the system, because the assumption (3.21) 
is satisfied here as shown in Appendix. As we calculate the one-loop contributions 
to two-point functions, we need Langevin equation to order a, which is 

K = d^a + a^ad^,^pd^,^p - a{N - l)(5"(0)$a + abN6''{0)^a + U 

Oi 

Ca = ((5a/3 - -$a$/3)r//3, (4.14) 

where 7]^ is white noise. We calculate two-point functions in one-loop order and 
show that there remains no (5*^(0) divergence. Integral equation corresponding to 
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eq.(4.14) is 



^a{x,t)^ j drG{x,t-T) r]a{x,T) - ^^a{x,T)^p{x,T)r]p{x,T) 

+a^a{x, T)d^^(3{x, T)d^^p{x, r) +a[bN-{N- l)^S\0)^a{x, r)] , (4.16) 
where 

j^-^ exp{-ikx)e{t)exp{-kh), (4.17) 
which can be solved by iteration as eq.(4.8). Eq.(4.16) leads to the following vertices 

Fig. 3 

In Fig.3 we denote G by a line, 770, a cross or an encircled cross and q;(5"^(0) a bullet, 
respectively. The above vertices contribute to the propagator shown diagramatically 
in Fig.4. 

Fig. 4 

In one-loop order the 5*^(0) -divergent contributions from Figs. (4c), (4d) and (4e) are 

7rT^^MS''(k + k') j <Pp + 2ml div. , 
(k^ + /c' )^ J 

With ^(0) = h, the sum of the above contributions is zero, which coincides with the 
result of the constrained system. It seems that in 0(N) non- linear a model both 
Langevin equations (2.19) and (3.1) lead to the correct results. 
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5. 5. Conclusion 



We constructed Langevin equation (2.19) for constrained system. In the deriva- 
tion of (2.19), the improved treatment of the stochastic consistency condition for 
constraints was essential. The Langevin equation can be applied to any system 
obeying (2.1)-type constraint. From the equation we directly derived Fokker-Planck 
equation and obtained eq.(2.22) as the equilibrium distribution, which coincides 
with the one obtained in path integral method. Owing to Ito calculus, the Langevin 
equation contains (5"'(0)-type singular terms and we showed explicitly, in 0{N) non- 
linear a model, that the singular terms are neccessary to the cancellation of (5"(0) 
divergences in one loop order. 

Furthermore, we ascertained that eq.(2.19) is equivalent to GCT-covariant and 

VR-invariant Langevin equation (3.1). We applied eq.(3.1) to 0(N) non-linear a 
model and showed that the singular terms in GCT-covariant and VR-invariant 
Langevin equation are also necessary to the cancellation of S^{0) divergences in 
one- loop order. 

6. Acknowledgments 

We thank Prof.S. Kawasaki for careful reading of this manuscript and also thank 
Dr.A.Nakamura for valuable discussion. 



16 



7. Appendix 



Here we show that assumption (3.21) is satisfied in 0(N) non-Hnear a model. We 
choose polar coordinates as new variables Q'^'s in section 3 and immediately recognize 
that e^jSI BG^i vanishes. If the assumption (3.21) is ascertained in the above special 
case, it is satisfied for any {Q^} because of the covariance of the assumption. 

We start with a set of variables {q'j}(i=l,2,- ■ -jN) obeying (qiqi)^ — rdvc = 
constant) and introduce a new set of variables {Q'^} (//=1,2,- • -jN). In accordance 
with section 3 

Q"" = Q"" = {q^ q^y^ - rc. 
If we choose N-1 angles of N dimensional polar coordinates as Q^'s {A=l,2,- • -jN-l) 

qi = rc cosQ^ 

q2 — rc sin Q cos Q , 

qN-i = rc sin sin • • • sin Q^~^ cos Q^~^ 
qN = rc sin Q sin Q • • • sin Q sin Q 



5 

N-1 



we obtain the vielbein e'i 



A dQ^ cos sin (5^^+^ • • • sin ^cosQ^ , , , 
oqi rc sm • • • sm ^ 

ei = 0, {A>i), 

. ^ - sin . ^ . 

' rc sinQi---sing^-i' ^ ^' 



e 



A 
N 



and e"- 



I 



cosQ^ sinQ^+i---sin(5^-i 
Vc sinQ^ • • • sinQ"^"^ 



= = sinQi • • • sing*-i cosQ\ {i N), 

eV = singi---sinQ^-\ 
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The vielbein satisfies (3.5) 



Vielbein leads to metric 



eYi = 0. 



gAB = ^A^B ^ I ^AB 

* * sin^ ■ ■ ■ sin^ Q^~^ 

By straightforward calculation we obtain 

eiVse^ = 0, 

which holds also for other choices of {Q^} due to the covariance of the equation. As 
we proved that assumption (3.21) is satisfied in 0(N) non- linear a model, we are 
allowed to use the equation (3.1) in section 4. 
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Figure Captions 

Fig.l Vertices from eq.(4.8) 

Fig.2 5"(0)-divergent diagrams contributing to two-point function 
Fig.3 Vertices from eq.(4.16) 

Fig. 4 Propagator including vertices shown in Fig.3 
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